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1. Introduction 

Spintronics as an emerging field of pliysics lias attracted considerable attention in recent 
years and has developed into various inter-related branches covered in this topical 
issue. Spintronics is devoted to employ the spin degree of freedom for information 
storage and as another means for extending the functionality of electronic systems. 
Semiconductor spintronics, as one subfield, is guided by the idea of combining concepts 
of spin electronics with the established techniques and advantages of semiconductor 
physics and nanostructures. Up to now, their properties used mainly rely on the 
charge degree of freedom alone. Many ideas for employing spin-polarized currents have 
been put forward since the seminal propsal by Datta and Das for a spin transistor [1] 
based on spin precession controlled by an external electric field through spin-orbit (SO) 
coupling [2]. These proposals usually require spin injection, more generally, the creation 
of spin-polarized particles in these materials. Spin injection from a ferromagnetic metal 
source into semiconductors is hindered by a fundamental obstacle originating from 
the conductivity mismatch between these materials [3]. Though this problem may be 
partially circumvented, for instance at low temperatures by tailoring dilute-magnetic- 
semiconductor/semiconductor interfaces [4, 5] enabling considerable spin-polarization 
ratios of the order of 90% [4], building all-semiconductor sources of spin-polarized 
electrons is still a challenge. 

Alternatively, several techniques to intrinsically create spin currents in non- 
magnetic systems have been put forward: Non-equilibrium spin-polarized currents have 
been created in two-dimensional electron gases (2DEGs) reahzed in zinc-blende-based 
heterostructures by means of optical pumping [6]. The irradiation of the 2DEG with 
circularly polarized light results in a spin photocurrent caused by the non-uniform 
distribution of the photoexcited carriers in k-space owing to optical selection rules and 
energy and momentum conservation. The recently proposed intrinsic spin-Hall effect 
in j9-doped [7] bulk systems and in 2DEGs [8] offers the principle possibility for spin 
current generation and manipulation in high mobility semiconducting systems. The 
combined effect of an applied electric field and the torque induced by SO coupling tilts 
the spins out of the 2DEG plane. This leads to spin accumulation at both sides of 
the sample in the direction perpendicular to the applied electric field. Spin-Hall effects 
have been observed experimentally both in semiconducting systems by optical detection 
techniques [9, 10, 11] and in metallic systems with an all-electrical setup [12, 13]. 

In the context of mesoscopic physics, further concepts such as adiabatic spin 
pumping [14, 15, 16] and coherent spin ratchets [17, 18] have been proposed for 
generating spin-polarized currents. This can be achieved by exploiting the magnetic 
properties of the semiconducting material, i.e. intrinsic spin-orbit interaction or the 
Zeeman coupling to external magnetic fields. 

Adiabatic quantum pumping involves the generation of a directed current in the 
absence of a bias voltage by periodic modulations of two or more system parameters, 
such as, e.g., the shape of the system or a magnetic field [19]. The spin analog of the 
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charge quantum pump can be achieved through an external magnetic field [15] or spin- 
orbit interaction [16] in order to filter the pumped current. For the case of an external 
magnetic field, it has been experimentally confirmed that under specific circumstances 
a spin current can be extracted from spin-dependent conductance fluctuations without 
accompanying net charge flow [20]. 

Ratchets [21] are generally systems with broken inversion (left/right) symmetry that 
generate directed net (particle) currents upon external AC driving in the absence of a 
net (time-averaged) bias potential. Thereby they have much in common with current 
rectifiers, though there are differences, in particular in the dissipative case [21]. The 
theoretical concept of ratchets, originally introduced for classical dynamics, was later 
extended to the quantum dissipative regime [22] . As a main feature, which distinguishes 
them also from rectfiers in the usual sense, quantum ratchets exhibit current reversal 
upon changing, e.g., the temperature or energy. Such quantum ratchets were 
experimentally realized in semiconductor heterostructures by demonstrating directed 
charge fiow in a chain of asymmetric ballistic electron cavities in the low-temperature 
regime, where the dynamics was close to coherent [23] . In this and further experimental 
investigations [24] also the charge current reversal phenomenon has been demonstrated. 

Very recently, we have proposed the generalization of the ratchet mechanism, 
extensively explored for particle motion, to generate directed spin currents. Spin 
ratchets require a coupling to the electron spin which, e.g., can be provided via spin- 
orbit interaction or external (non-uniform) magnetic fields. In Ref. [17] spin-orbit 
ratchets have been considered in the coherent regime. There a proof of principle for 
a net ratchet spin current (in absence of an average charge current) has been given 
and confirmed by numerical calculations for experimentally accessible parameters for 
GaAs-based heterostructures. 

In the present paper we further investigate the possibility of employing Zeeman 
ratchets for spin current generation, i.e. by considering mesoscopic conductors with 
a spatially varying Zeeman term and subject to an AC bias. In the simplest case 
of one-dimensional motion, and assuming preserved spin states, spin-up and -down 
electrons will experience a Zeeman term, ±{g* /2)fiBB{x), where B{x) is the external 
non-uniform magnetic field, g* the effective gyroscopic factor and /xb the Bohr magneton. 
li B{x) 7^ B[—x), then the different spin species experience opposite asymmetric Zeeman 
ratchet potentials. In close analogy to the particle ratchet mechanism described above, 
the different spins are expected to be predominantly driven into opposite directions 
upon external driving, resulting in a spin-polarized current. This mechanism has been 
studied and confirmed in Ref. [18]. 

In the present work we relax both assumptions of one-dimensional motion [25] 
and conserved spin directions, e.g. the notion of two independent spin species, and 
explicitly include spin-fiip effects. Moreover, we work out how such spin ffip processes 
can be invoked to engineer and tune ratchet spin conductances. To this end we consider 
spin ratchet effects of a system consisting of a two-dimensional quantum wire embedded 
in 2DEG and subject to an AC bias in between two ohmic contacts. The non-uniform 
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Figure 1. Ferromagnetic stripes (magnetization direction given by red arrows) on 
top of a semiconductor heterostructure that harbors a two-dimensional electron gas 
(indicated by white lines) containing a quantum wire (black) of width W. a) Setup A: 
two stripes with antiparallel in-plane magnetization, see Sec. 3; b) setup B: one stripe 
with out-of-planc magnetization, see Sec. 4. 



S-field is created by the magnetic fringe fields of ferromagnetic stripes patterned on the 
semiconductor heterostructure [26]. 

This paper is organized as follows: In Sec. 2 we introduce the model for the spin 
ratchet. We thereby specify the driving of the ratchet and the evaluation of the net 
charge and spin currents. In Sec. 3 we study transport through the quantum wire subject 
to the magnetic fringe fields of two ferromagnetic stripes with antiparallel magnetization 
perpendicular to the quantum wire in the plane of the 2DEG. This configuration (setup 
A in Fig. la) allows us to study the transition from decoupled to strongly coupled spin 
states and its implications on transport and thus also on the ratchet currents. In Sec. 4 
we then investigate the conductor subject to the fringe field of a single ferromagnetic 
stripe magnetized perpendicular to the plane of the 2DEG (setup B in Fig. lb)). Using 
symmetry arguments and numerical calculations we demonstrate that the two setups 
introduced in Sec. 3 and 4 act as spin ratchets. After summarizing in Sec. 5, we close the 
paper with an appendix on the general derivation of an expression for the spin current 
within the framework of the multi-terminal Landauer-Biittiker formalism and a second 
appendix including the derivation of symmetry relations for the transport properties at 
finite applied bias. 

2. Model and formalism 

We consider a quantum wire in the x-direction embedded in a 2DEG in the (x, y) plane. 
The system is subject to a non-uniform magnetic field B{x,y), due to the fringe fields 
of ferromagnetic stripes patterned on top of the 2DEG (see Fig. 1). Their deposition on 
a semiconductor heterostructure can be accomplished with electron beam lithography 
and lift-off techniques [27, 28]. Near-surface 2DEGs can be fabricated to lie only a few 
tens of nanometers beneath the surface [28], thereby achieving magnetic field values of 
up to 0.5T with thin ferromagnetic films [29]. The magnetic fringe field of a ferromagnet 
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with homogeneous magnetization M is given by [30] 

B(r) = -tlv i ia' (2.1) 



where the integration runs over the surface S of the ferromagnetic stripe, and u{f') 
is the unit vector normal to the surface of the stripe at position r Accordingly, the 
corresponding vector potential in the Coulomb gauge, V • >1 = 0, is given by [30] 

r/^N Mo / , ,Mxu(f') 

We model the wire, where electron transport is assumed to be phase coherent, by the 
single-particle Hamiltonian 

^0 = ^ ^-B{x,y) ■ a + V{y), (2.3) 

where g* is the effective gyroscopic factor, m* the effective electron mass, /ib the Bohr 
magneton and a the vector of the Pauli matrices. The term V{y) denotes the lateral 
confining potential defining the quantum wire. Orbital effects due to the magnetic 
field are accounted for by the vector potential A{x,y) in Il{x,y) = p — eA{x,y). 

— * 

Spin effects in transport through the wire enter via the Zeeman term {g*ii^/2)B ■ a 
coupling the spin degree of freedom to the external magnetic field. For a proper 
treatment of the spin evolution, the inclusion of the full magnetic field profile is 
mandatory [31]; disregarding [32] one of the magnetic field components Bi may lead 
to an oversimplification of the problem. 

To obtain a significant spin ratchet effect materials with a large g* factor are most 
suitable. In this respect dilute magnetic semiconductors (DMS) represent a promising 
class of materials. Recent measurements have shown values of g* > 100 [33] where 
in addition to the intrinsic g^-factor, an additional contribution to g* appears, owing 
to exchange coupling among the electron spins and the magnetic ions present in the 
DMS [34]. These materials with a large g* factor can also exhibit large SO couphng 
values. The working principle of a spin ratchet based on SO coupling has been already 
investigated in Ref. [17]. As one result of this study a spin-ratchet effect can only occur 
if there exists the possibility to mix different transverse channels of the wire. This can 
only happen when more than one open channel is taking part in the transport and 
electrostatic barriers induce this mixing of different bands. As we will show below, the 
spin-ratchet effect created by the setup presented in this paper does not rely on this 
condition and is already present when only the first conducting channel is opened. Due 
to this fact and the absence of any electrostatic barriers responsible for the mixing of 
different subbands we disregard here the SO couphng and focus on the effects due to 
the presence of the magnetic stripes. 

In the present work where we consider disorder-free ballistic motion, we refer to 
nonmagnetic high-mobility semiconductors. To be definite we chose throughout the 
paper parameters for InAs 2DEGs with typical values of m* = 0.024mo[35], where mg 
is the free electron mass, and \g*\ = 15. InAs is well suited due to its large g* factor and 
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the property that InAs 2DEGs can be fabricated close to the surface where the magnetic 
stripes are located. We assume that the stripes possess a magnetization /XqM — 3T. 

The charge current Iq through the wire is evaluated within the Landauer approach. 
For coherent transport in a two-terminal device the current can be expressed as 

e f°° 

Ic = dE [fiE; /xl) - f{E; /.r)] TiE) (2.4) 

in terms of the quantum probability T{E) for electrons with energy E to be transmitted 
from the lead at higher to the lead at lower potential. In Eq. (2.4), /(i5;/XL/R) is the 
Fermi function for the left /right lead with chemical potential /xl/r- 

The spin current Isix) passing a cross section {x = const) is given, for a 
wavefunction ^(x,y), by 

/s(a;) = J dy^*{x,y)Js'if{x,y). 

Here we use the most common definition [36] of the spin current operator Js which, 
with respect to an arbitrary quantization axis u, reads 

inside the leads. The partial derivatives in (2.5) act on expressions to their right and 
left (indicated by the arrows) . 

Contrary to the charge current that obeys a continuity equation, the spin current 
can take different values if evaluated in the left or in the right lead. This usually happens 
in systems where the Hamiltonian does not commute with the Pauli matrices a giving 
rise to a torque inside the scattering region [37], which can change the spin state of the 
electrons. For this reason we will explicitly label the lead, where we compute the spin 
current. Although there is some freedom in the choice of the spin current operator [38], 
here we evaluate the spin current inside the magnetic field free leads, where a tourque 
term in the continuity equation for spin is absent. Therefore the spin current inside 
the leads is a well defined quantity, which can be measured in principle. As derived 
in Appendix A (Eq. A. 9) the corresponding spin current in the right lead reads 

1 r°° 

7s = dE [f{E-t,^) - f{E;i,^)]Ts{E) , (2.6) 
with the spin transmission probability defined as 

Ts{E)^ J2 [T^AE)-T.AE)]- (2.7) 

a=±l 

Here T^'^a is the probability for an electron with initial spin state a to be transmitted 
from the left lead into the spin state a' inside the right lead, see Eq. (A. 7) in Appendix 
A. To obtain the corresponding spin current in the left lead one has to replace T^r^^{E) 
in Eq. (2.7) by the corresponding probabilities T^, ^{E) for transmission from the right 
to the left lead. 
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Ordinary particle ratchets give rise to a net drift motion of particles in one 
preferential direction upon AC driving without net bias (rocking ratchet). Below we 
will generalize this concept to induce spin-dependent ratchet currents correspondingly. 
The AC driving can be considered as adiabatic, since the timescales for the variation of 
an external bias potential are long compared to the relevant time scales for charge 
transmission through the device. For a proof of principle we assume an adiabatic 
unbiased square- wave driving with period to. The system is periodically switched 
between two rocking conditions, labeled by bias ±Uo {Uq > 0). The electro-chemical 
potential //l/r of the left /right reservoir is changed periodically in time according to 

f ep ± Uo/2 for 0<t < to/2 

[eF=F^o/2 for to/2 < t < to 

^•6. /iL/R(t) = AtL/R(t + to)- In adiabatic limit considered, the system is assumed to 
be in a steady state between the switching events. Then the ratchet charge and spin 
currents inside the wire are obtained upon averaging Eqs. (2.4,2.6) between the two 
rocking situations 

(7c(£f, Uo)) = \ [/c(£f, +^o) + /c(£f, -^o)] (2.9a) 



= - — dE Af{E; £F, Uo)AT{E; Uo) , 



(/s(£f, Uo)) = \ [/s(£f, +^o) + /s(£f, -Uo)] (2.96) 



1 

8"7r Jo 



dE Af{E;ep,Uo)ATs{E;Uo), 



where 



Af{E- ep, Uo) = f{E- £f + Uo/2) - f{E- - Uo/2), 

AT{E;Uo) =T{E;+Uo)-T{E;-Uo), (2.10) 

ATs{E; Uo) = Ts{E; +Uo) - Ts{E; -Uo). 

An extension to an adiabatic harmonic driving is straight forward. 

In linear response, the ratchet spin current (Is), Eq. (2.96), vanishes because 
ATs{E; Uq = 0) = 0, see Eq. (2.10). Hence we must consider the nonlinear regime to 
obtain a finite net spin current. Since we consider nonlinear transport ignoring inelastic 
processes, we can write the currents Eq. (2.4) and (2.6) as energy integrals over the 
transmission. To model a finite voltage drop across the two leads, we add the term 

nu = Ug{x,y;U) (2.11) 

to the Hamiltonian (2.3), where for the square wave driving considered here, U takes the 
values ±f/o respectively. Here the function g{x,y; U) describes the spatial distribution 
of the electrostatic potential inside the mesoscopic system and is generally obtained 
through a self-consistent solution of the many-particle Schrodinger equation and the 
Poisson equation [39, 40]. However, here we employ heuristic models for g{x,y;U), 
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assuming that the vohage primarily drops in regions, where the magnetic field strongly 
varies spatially. This model is based on the fact that the corresponding Zeeman term 
acts as an effective potential barrier, and takes into account that a more rapid potential 
variation leads to enhanced wave refiection and hence to a steeper local voltage drop [41] 
(details will be given in the following sections). 

In order to numerically evaluate the transport properties of the system, the stationary 
Schrodinger equation {H — E)^{r) = is discretized on a square lattice, yielding a 
tight-binding representation of H. This is then used to calculate the elements of the 
scattering matrix of the system via lattice Greens functions and a recursive Greens 
function algorithm [42]. 

3. Setup A: Two ferromagnetic stripes with longitudinal magnitization 

In this section we investigate the spin-dependent transport properties (and thereby also 
the operabihty as a ratchet) of a quantum wire in the plane of the 2DEG subject to 
the magnetic field of two stripes with opposite, longitudinal in-plane magnetizations, 
M — ±My, arranged perpendicular to the wire. This setup A is shown in Fig. la). 
The Hamiltonian of the system is given by Eq. (2.3) with Ay — (in the Coulomb 
gauge), i.e. Uy = py. For the following analysis we chose a confinement potential V{y) 
such that the wire of width W is displaced by a shift yo with respect to the symmetric 
configuration, see Fig. 1. 

For sufficiently narrow wires {W <^ b) and small displacement {yo -C b), the energy 

A — * — * 

scales of the magnetic field contributions to Ho in Eq. (2.3) containing B and A are much 
smaller than differences between energy levels E^n — of different transversal modes 
|m) and \n). Therefore transitions between different transversal modes are strongly 
suppressed, and we consider the case of only one open mode. Higher modes \n) mimic 
the behavior of the first mode up to an energy offset En — Ei. 

Evaluating Eq. (2.1) for the combined magnetic field of two stripes centered around 
(x = 0,1/ = 0), setup A, we obtain the following symmetry properties for the 5-field 
components of this configuration: 

B^{x,-y,z) ^ -B^{x,y,z), (3.1a) 

By(x,-y,z) ^ By(x,y,z), (3.1b) 

Bz{x,-y,z) = -B^{x,y,z) . (3.1c) 

In particular, Eqs. (3.1a) and (3.1c) imply vanishing magnetic field components B^ and 
Bz for y = 0. Therefore, we use y as the spin quantization axis for the considerations 
below. Those symmetries also have an important implication for the spin dynamics of 
the system. For a confinement potential that is symmetric upon refiection at the {x, z)- 
planc, V{—y) = V{y), as realized for a symmetric confinement centered around yo = 0, 
spin-up and spin-down eigenstates within the same transversal mode n are decoupled. 
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Figure 2. Magnetic field components (dasli-dotted green line), By (solid red), —By 
(dotted blue) and B^ (dashed black) in the plane of the 2DEG at fixed y = 200nm for 
setup A, Fig. 1(a) and parameters given in the text. 




The relevant matrix element {n,a\ a ■ B{x,y) \n, 
vanishes, 



dy\Xn{y)f B,{x,y) -i(rB^{x,y) 



-a), responsible for the spin mixing, 



0, 



(3.2) 



since the integrand is an odd function of y due to Eqs. (3.1a,3.1c) and the fact that 
transversal modes obey Xn{—y) = (— l)"~^Xn(l/)- However, for finite values of yo, the 
coupling (3.2) does not vanish anymore, and spin flips can arise with a significant effect 
on electron transport, as we will demonstrate in the next paragraph. 



3.1. DC transport 

Before evaluating the average charge current (2.9a) and spin ratchet current {2.9b) it is 
instructive to analyze the DC transport properties of setup A. To this end we chose as 
realistic parameters for the geometry (see Fig. 1(a)) W = 120 nm, a = 600 nm, b = 2/im 
(thus W/b = 0.06), c = 200 nm, d = 600 nm and e = 100 nm. For this parameter set and 
I/O -C 6 the magnetic field component By{x,y) is approximately constant in ^/-direction, 
i.e. By{x,y) ~ By{x). It possesses a much larger maximum value than the other 
components B^ and B^. In Fig. 2 we show the x-dependence of the overall magnetic 
fringe field of the two-stripe setup A for fixed y = 200nm. 

Figure 3 shows the total transmission T{E) in linear response [Uq —>■ 0) for 
energies within the first transversal subband for different values of y^. For yo = 
the spin eigenstates decouple owing to Eq. (3.2). The energy where the first 
transversal mode opens is shifted to ~ i?i + [/b due to the Zeeman barrier of height 
Ub = g*fJ^B'^8.x[By{x,y = 0)]/2 present in the wire. However, for increasing yo when 
spin fiips can take place, an additional plateau builds up at energies Ei < E < Ei + Ub 
approaching T ^ 1. In the inset of Fig. 3 the spin- resolved transmission probabilities 
To-',o- are depicted for yo = 100 nm. We identify T^-{E) as the sole contribution to the 
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Figure 3. Total transmission T{E) in linear response for values of yo = nm to 
200 nm in steps of 50 nm from bottom (black line) to top (brown line). The dotted 
black line indicates the transmission for B — 0. Inset: spin-resolved transmission 
probabilities T^'^a, Eq- (A. 7), for yo = 100 nm. 



total transmission T{E) for Ei < E < Ei + Ub- Thus the appearance of the additional 
plateau is a consequence of the mixing of the spin states. For energies well above the 
barrier both, spin-up and spin-down electrons are fully transmitted. 

The main features in the numerically calculated transmission in Fig. 3 can be 
understood using a heuristic model. It is based on the fact that in the region close to 
X = (see Fig. 2) spin flips predominantly take place, since By{x,y) = vanishes at 
X = 0, and spin- up and -down states of the same transversal mode are nearly energy 
degenerate. 

In the following we consider stepwise (positions labeled in Fig. 2) the spin evolution 
along the wire for unpolarized electrons entering the system with energy Ei < E < 
El + Ub: 



Unpolarized electrons (equal number of spin-up and spin-down particles) are 
injected from the left reservoir. 

Spin-up electrons are completely reflected at the Zeeman barrier (indicated by the 
solid red line), while spin-down electrons experience a potential valley (blue dotted 
line) and can pass. 

A fraction of the spins is flipped from down to up due to a finite Bz{x) close to 
X = 0. 

(4) Spin-down electrons are completely reflected while spin-up electrons pass. 

(5) Only spin-down electrons from the left lead reach the right lead, after undergoing 
a spin flip. 

Hence this mechanism leads to T+ + = T_ + = T_ _ = and „ 7^ 0. Although 
the model can explain the basic features of the transmission curves shown in Fig. 3 
fairly well, it cannot account for the details in the functional dependence of T[E) for 
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energies below the barrier which reflects further quantum effects present, e.g. resonant 
tunnehng processes. An analysis of the spin-resolved transmission probabilities in the 
opposite rocking situation, /^l < /^r, where electrons flow from right to left, shows 
that transmitted particles are oppositely spin-polarized compared with the former case. 
Correspondingly, T!_ _,_ is the only non-zero component of the spin- resolved transmission 
for energies E < E < Ei + U-q- 

The above analysis demonstrates that the magnetic field components perpendicular to 
the dominant one, even if they are small, can significantly alter the transport properties 
of the system. In the present case, disregarding and would have resulted in a 
vanishing transmission for E < Ei + Ub- 

3.2. AC transport 

We now investigate the rectification properties of setup A upon applying the AC driving 
given in Eq. (2.8). We first specify how we obtain the drop of the electrostatic potential 
g{x,y; U) across the system for a finite applied bias. Based on a heuristic model used 
in Ref. [24] we assume that {dg{x) / dx) oc \{d/dx)By{x,y — 0)\ in the central scattering 
region, —L/2 < x < L/2, yielding 

1 f':^,^dx\{d/dx)By{x,y^O)\ 

gix) = ptt; , (3.3) 

2 J^i%dx\{d/dx)By{x,y^O)\ 

while we fix g{x) to ±1/2 inside the left (right) lead. For the full system Hamiltonian 
H — Ho + Hu: Eqs. (2.3,2.11), at finite bias Uq we evaluate the expressions (2.9a,2.9&) 
for the average charge and spin currents. If orbital effects due to the perpendicular 
magnetic field Bz are negligible, i.e. Ay ~ 0, the total Hamiltonia Til is invariant under 
the symmetry operation V = R^Ruf^z- Then, as shown in Appendix B, the relation 
(B.4) between S'-matrix elements holds true. Squaring S'-matrix elements in Eq. (B.4) 
and summing over channels (ncr) G L and (nV) G R yields T{E-^±.Uq) = T'{E;^Uo). 
This relation, together with the relation T'(E; tUq) = T{E] ^Uq) due to unitarity of 
the (S'-matrix, leads to (/c) = 0. A vanishing average charge current is in line with 
symmetry considerations for charge ratchets and coincides with a numerical analysis for 
the parameters used here. 

However, on the other hand, the symmetry considerations imply that the average 
spin current can take finite values. To conflrm this numerically and to get an idea 
of its magnitude we calculate the ratchet spin current in the right lead according 
to Eq. (2.96). Figure 4 shows the differences in spin transmissions, ATs{E), for the 
two rocking situations as a function of energy for moderate flnite applied bias voltage 
Uq — O.IC/b- For yo — 0, where transitions between spin-up and spin-down states of the 
same transversal subband are absent, the system is comparable to the devices studied 
in Ref. [18]. There it was shown, in analogy to the case of charge rectification [43], that 
for conserved spin eigenstates the ratchet effect stems from different maximum values 
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Figure 4. Ratchet spin transmission ATs for spin quantization axes x (Panel a), 
y (Panel b) and z (Panel c) as a function of the injection energy for displacement 
Uo ^ nm (solid black line), 100 nm (dashed blue) and 200 nm (dash-dotted red) for 
bias potential Uq = O.IC/b and g{x) specified in Eq. (3.3). For comparison, ATs at 
2/0 = nm (brown triangles) and j/o = 100 nm (brown circles) is shown for a linear 
voltage drop g{x) — —x/L across the central scattering region with bias potential 
Uo ~ 2.3C/0. 



Amax(t^o) of the effective potential landscape in the two rocking situations, 

Amax(f/o) = m&xpcsA^^ +^o)] " max[Ues,a{x, -Uo)], (3.4) 

with 

f/eff,<,(x, ±Uo) = ±Uogix) + {a/2)fij,g*By{x, y = 0) . (3.5) 

This mechanism explains the rectification related to ATs-y{E;yo = 0) (solid black line 
in Fig. 4(b)) and its functional dependence for yo = 0. There spin flips are absent, and 
hence T+ _ = T_ + = 0. For increasing yo the magnitude of ATs-^{E) decreases, while, 
at the same time, ATs;x{E) and ATs-^iE) grow and take finite values. Thus we can 
summarize that the ratchet effect survives in the presence of mixing of different spin 
states. However, as apparent from Fig. 4, the vector of spin polarization is no longer 
aligned along y for finite yo. 

We further study how sensitively the observed effect depends on the particular form 
of the voltage drop g{x). To this end in Fig. 4 we additionally show ATs for a linear 
voltage drop model g{x) = —x/L (brown circles/triangles) inside the central region 
{—L/2<x<L/2), where the bias voltage IIq was chosen such that the maximum value 
of the respective effective potential (3.5) was the same for both voltage drop models: 
max[[/eff,^(a;, ±[/o)] = max[f/efr,a(a;, ±f^o)]- Comparing ATs{E;yo) for g{x) at |/o = 
(brown triangles) and 100 nm (brown circles) with the respective function ATs{E; y^) for 
g{x) we observe no significant difference in their functional dependence on although 
the magnitude of the overall bias is different for both models {Uq ~ 2.?>Uq for the 
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Figure 5. Bias voltage dependence of the ratchet spin conductance {Is;u) {e / Uq) at 
zero temperature, k^T = 0, for a Fermi energy ep = Ei + Ub- Results are shown for 
tjq = and polarization axis u = y (black solid line), and for yo — 100 nm and u ~ x 
(red dashed line), y (green dash-dotted hne) and z (blue dotted hne). 



curves presented in Fig. 4). Therefore we can state that ATs{E] yo) rather depends on 
the difference (3.4) in the maximum values of the effective potential than the actual 
distribution of the electrostatic potential in the mesoscopic conductor. A study of ATs 
in the other lead shows very similar results to the ones presented in Fig. 4. 

In Fig. 5 we finally display the ratchet spin conductance (/s(£:f, f^o))(e/f/o)) 
Eq. (2.96), which shows a nearly linear dependence on the bias voltage Uq, i.e. 
ATg(£^;i/o) oc Uq. This is in line with the above analysis showing that ATg(i?;?/o) 

'—'max 

[Uq) and Amax(t^o) oc Uo. 

4. Setup B: a single ferromagnetic stripe with transverse magnetization 

In the following we investigate the possibility of generating a spin ratchet effect 
using the magnetic field profile of a single stripe magnetized in the (x,z)-plane as 
shown in Fig. 1(b). For b ^ W the evaluation of Eqs. (2.1,2.2) for a magnetization 
M = M^x + MJ yields 





/ 




[x,z) 


\ 




/ 





\ 


B{r} = 











, A{r) = 




Ay(^X, Z) 






v 




{x,z) 


) 




\ 





/ 



The Hamiltonian of this system then reads 

+ Viy). 

The 2DEG is located e = 100 nm below the surface of the semiconductor heterostructure 
(see Fig 1(b)). The extension of the stripe in x-direction is chosen to be a = 600 nm, 
infinite in y-direction and c = 200 nm in z-direction. For the analysis below we chose 
a stripe magnetization M = Mz. The corresponding magnetic field in the plane of the 
2DEG is depicted in Fig. 6. Results comparable to those presented below are obtained 
for a stripe magnetized in x-direction. 



Ho 



Px + . 

2m,* 



.X 



Bx{x)(T^ + B^{x)a^ 
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X 



Figure 6. Magnetic field components Bx{x) (solid black line) and Bz{x) (dashed 
red line) in the plane of the 2DEG produced by a ferromagnetic stripe (a = 600 nm, 
b oo, c = 200 nm, e = 100 nm, see Fig. 1(b)) with magnetization M = Mz. 




Figure 7. Total transmission T{E) (solid black line) and absolute value of the spin 
transmission |rg| (dashed red line) as a function of energy for a wire underneath a 
single ferromagnetic stripe, see Fig. lb and text. For comparison, the dotted brown 
staircase function shows the transmission in the absence of a magnetic field. 

DC transport 

In Fig. 7 we show the total transmission T{E) for a quantum wire of width W = 120nm 
subject to the 5-field, Fig. 6, in hnear response {Uq — > 0). In addition to the steps 
at even values of T{E) ^ 2,4,6... due to the successive opening of the transversal 
modes at energies = [h'^Tr'^/{2m*W^)]n'^, additional plateaus appear at odd values 
of T{E) ^ 1,3,5... close to the energies En- As in Section 3, we can attribute 
these features to the lifted spin degeneracy due to the Zeeman field, since also the 
width of these plateaus corresponds to twice the absolute height of the Zeeman barrier 
Ub = (fi'*/2)/iB niax[ ] inside the wire. 

In Fig. 7 we furthermore plot the absolute value of the spin transmission, 
|Ts| = a/(7s;xO^ + (^S;?/)^ + (^S;^)^) which approachcs unity at energies of the additional 
plateaus. A closer look at the spin- and mode-resolved transmission probabilities reveals 
that the transmission of the highest occupied transversal subband is completely spin 
polarized at the plateaus, while the lower modes are fully transmitting spin-up and 
spin-down particles. Similar results were reported in Ref. [31] for a stripe magnetized 
in the x-direction. 
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Apart from the spin effects due to tlie Zeeman term, the vector potential component 
Ay, affecting the orbital dynamics of the electrons due to the perpendicular magnetic 
field Bz, influences the electron transport. In a classical picture, Bz forces electrons to 
move on segments of cyclotron orbits in the plane of the 2DEG. Therefore, the kinetic 
energy in the direction of motion is reduced resulting in a shift of the energies where 
the transversal modes open towards higher values [44]. This is visible in Fig. 7, when 
comparing the total transmission with (solid black line) and without (dotted brown line) 
magnetic field. 

4-2. AC transport 

As for setup A we employ a heuristic model for the voltage drop inside the mesoscopic 

— * 

conductor, assuming dg{x)/dx oc \{d/dx)\B{x)\\ inside the central region {—L/2<x< 
L/2) yielding 



I-L/2 


{d/dx)\B{x)\ 




{d/dx)\B{x)\ 



However, before we numerically investigate the AC ratchet transport properties, we 
exploit certain symmetries present in the system to simplify the expressions for the 
average net charge (2.9a) and spin currents {2.9b). For the magnetic field profile 
produced by a stripe magnetized in 2;-direction it is straightforward to show from 
Eqs. (2.1,2.2) that the following symmetry relations hold true (see also Fig. 6): 

B,{-x) = -B,,{x) , Bzi-x) = Bzix) , 

Ai-x) = -Ai^) ' 9{-x) = -g{x) . 

Thus the Hamiltonian H — + Hu is invariant under the action of the operator 
V — —iCRxRu^^z-, where C is the operator of complex conjugation, inverses the x- 
coordinate, Ru changes the sign of the applied voltage (±f/o ^ T-Uq) and o"^ is the Pauli 
spin operator. Due to \H, P] = we are able to interrelate the transmission probabilities 
for both rocking situations as shown in Appendix B. Taking the square of Eq. (B.3) and 
summing over the transversal modes n G L and n' e R we obtain the following relations 
between the spin-resolved transmission probabilities in the two rocking situations: 

Here the superscript labels the angles of the spin quantization axis on the Bloch sphere 
(see Appendix B). Thus the ratchet charge current (2.9a) vanishes, 

(7c(£f,C/o)) = 0, 

and we can express the ratchet spin current (2.9&) through the transmission probabilities 
of a single rocking situation {e.g. -\-Uq): 

1 f°° 

fe(^F,t/o)) dEAf{E;eF,Uo)x 
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Figure 8. Averaged net spin transmission ATs{E;Uo) for spin polarization axis x 
(solid black line), y (dash-dotted red line) and z (dashed blue line) as a function of 
energy close to the lowest three transversal energy levels En {n — 1, 2, 3) for an applied 
bias voltage of Uq = O.IUb- 



X 



{hiy/zi^F, Uo)) 



T+,4E,+Uo)-T^AE,- 
dE A/(E;£p,f/o) 



-Uo) 



X 



1 

47r 



X 



iE,+Uo)-T^AE,+Uo) 



Figure 8 shows the ratchet spin transmission ATs{E; Uq) at a finite apphed voltage 
Uq = O.IUb for a wire of width W = 120 nm. This quantity is finite for energies where 
the DC transmission is spin polarized (see Fig. 7). Furthermore, the spin polarization of 
the ratchet spin transmission depends on the injection energy. This opens the possibility 
to tune it upon varying the Fermi energy. Comparing the ratchet spin transmission for 
n = l (Panel a) and n = 2, 3 (Panels b,c) we observe that its magnitude is significantly 
lower in the case where more than one transversal mode is conducting. This behavior 



is due to the mixing of different transversal subbands due to Ay 
effect we introduce 



To quantify this 



AT. 



S,max 



n] 



max 

-£'6[-Bn,-En + l] 



as a measure for the rectification in each single transversal mode. Figure 9 shows that 
ATs.max('^ = 2)/ATs,max('^ = 1) ^ 1 for cascs where the mixing due to Ay{x) is small, 
i.e. for a narrow wire and/or small magnetic field. However it decreases upon increasing 
fioM and/or W. Note that for setup A in the previous section the magnetic field inside 
the quantum wire was one order of magnitude smaller than for setup B here, thus 
yielding a comparable value ATs,max('^) for all subbands n. 

As for setup A, the ratchet spin conductance shown in Fig. 10 exhibits a linear 
dependence on the applied voltage. Thus we presume that the rectification mechanism 
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Figure 9. Ratio ATg inaxl"- = 2)/ATs niaxl"- — 1) as a function of (a) the 
magnetization of the ferromagnetic stripe and (b) the wire width. 



0.03 r 

K 0.02- 




Figure 10. Ratchet spin conductance (/s.«)(e/C/o) as a function of bias voltage Uq at 
k^T — and ep = E3 + 3Ub for the spin quantization directions x (black solid line), 
y (red dashed line) and z (green dash-dotted line). 

is the same as in setup A, although the spin dynamics is much more intricate. 
5. Conclusions and Outlook 

In the present work we have shown that the coupling of the electron spin to the magnetic 
fringe fields of ferromagnetic stripes via the Zeeman interaction can be used to generate 
a spin ratchet effect in a coherent mesoscopic conductor subjected to an adiabatic AC 
driving with finite bias. The proposed devices exhibit the appeahng property of creating 
a directed net spin current in the absence of an accompanying net charge transport. This 
key result has been demonstrated in numerical approaches for setups A and B, in the 
case of setup B, also analytically based on symmetry properties of the system. 

The generated spin current may be regarded as resulting from a rectification effect, 
however in a generalized sense: The direct analogue of a charge current rectifier would 
be a system generating a directed spin current out of a conductor with alternating spin- 
chemical potentials in the left and right reservoir. Spin ratchets, as the ones considered 
here, act differently as they convert an AC electrical bias into a net spin current. 

/^From our analysis we have identified the difference in the maximum values of 
the effective Zeeman potentials in the respective rocking situations as responsible for 
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the creation of the spin current. It has been shown that this rectification effect is 
almost independent of the actual distribution of the electrostatic potential in the biased 
conductor. Furthermore, the fact that, for the systems considered here, the difference in 
the maximum values of the Zeeman potential is crucial for the spin currents, implies that 
the magnitude of the spin current cannot be systematically increased upon increasing 
the number of magnetic stripes, e.g. in a periodic arrangement of stripes. We have 
checked this also numerically by adding an increasing number of stripes. 

In the preceding sections we presented results, when evaluating the spin current 
inside the right lead. However, as we noted both systems considered, setup A and B, 
are characterized by interesting symmetry properties. Those can now be used to directly 
extract the respective currents inside the left lead. If for setup A the component 
can be neglected, as it is appropriate for the parameters used in Sec. 3, the combined 
Hamiltonian of the system and the leads has been proven to be invariant under the action 
of the symmetry operation R^RyRuaz, while for setup B the Hamiltonian is invariant 
under the action of RxRuO'z- These symmetry properties are reflected in Eqs. (B.4) 
and (B.5) respectively, see Appendix B. Both relations lead to the general relation 

T:J^^Hi?,Tf/o)=T!;f)(E,±f/o) 

between the transmission probabilities in the two rocking situations. This relation allows 
for the following interpretation. If the transmitted electrons arc spin-polarized in one of 
the two rocking situations with direction of the polarization vector given by the angles 
{9, (f)) on the Bloch sphere, then in the other rocking situation the spin polarization 
vector of the output current (in the other lead) will be rotated around the z-axis by 
TT and thus points to (^,0 + tt). This property is a direct consequence of the lacking 
conservation of the spin eigenstates inside the wire. 

The physics of semiconducting materials characterized by a large g* factor is dominated 
by the presence of magnetic impurities, e.g. diluted magnetic semiconductors [34]. 
Therefore, in order to exploit the expected stronger rectification effect, we have to 
take elastic scattering off impurities into account. In particular, we plan to study how 
additional disorder alters the spin ratchet effects. 

Finally, since the heuristic model used for the distribution of the electrostatic 
potential in the conductor is convenient but not fully satisfactory, it would be desirable 
to calculate the charge density and the respective electrostatic potential inside the wire 
self-consistently. Work in this direction is in progress. 
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Appendix A. Calculation of the spin current in the Landauer-Biittiker 
formalism 

In this appendix we present a derivation of the expressions for spin current in the leads 
of a multiterminal coherent conductor within the framework of the Landauer-Biittiker 
theory [45]. To this end we consider N non-ferromagnetic contacts, injecting spin- 
unpolarized current into the leads. For convenience we use a local coordinate system for 
the lead under investigation, where x is the coordinate along the lead in the direction 
of charge propagation due to an applied bias in linear response and y is the transverse 
coordinate. Then the eigenfunctions inside a lead are given by 



(A.l) 



where the Xn{y) are the transverse eigenfunctions of the lead with the transversal 
eigenenergy En and S(cr) is the spin eigenfunction. The superscript ± of $ refers to the 
direction of motion in ±a;-direction with the wave- vector kx = \/2m*{E — En)/ h. For 
the derivation we use the scattering approach, where the amplitudes of the states inside 
the leads are related via the scattering matrix S(£'), determined by the Hamiltonian of 
the coherent conductor. Inside lead q a given scattering state 

^%{^,y) =Y.i<Am%nA^.y) + KAmE,nA^.y)) , 

((7 = ±), consists of incoming states entering the coherent conductor from contact q 
and outgoing states $~ leaving the coherent conductor into contact q. The amplitudes 
of incoming a^^ and outgoing waves 6^,^, are related via the equation 



N 



(A.2) 



j=l nGj (T=±l 

where the scattering matrix S{E) has the following structure for an N terminal system: 



S{E)^ 



[ v^'^E) 
t''^\E) 



t^'^iE) 
(E) 



K.2,2 



tl'^(E) \ 

t2'^(£;) 



(E) J 



\ t^'i(^) t^'^{E) ■ 

Here the sub- matrix r^'^[E) is a square matrix of dimensionality M^[E), correspond- 
ing to the number of open channels at energy E in lead j (already including the spin 
degree of freedom), which is connected to a reservoir with chemical potential /jj. The 
matrix r^'^[E) contains the scattering amplitudes of incoming channels of lead j being 
reflected back into outgoing channels of the same lead. The sub- matrix t*'^(£') is a 
M^{E) X M^{E) matrix that contains the scattering amplitudes for transmission be- 
tween incoming channels from lead j and outgoing channels of lead i. 
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The wave function of the scattering state inside lead i, where only the incoming 
channel (na) e j is populated (a^/^/ = Sj/jSn',n^a',(7): reads for j — i: 

[n'(j')& 

and, correspondingly, for j ^ i 

For a wave function ^{x,y) the spin current I^{x) passing a cross section {x = const) 
of a lead is given by: 



(A.5) 



Here we use the most common definition of the spin current operator [36], which with 
respect to an arbitrary quantization axis u takes the following form inside the leads: 



7 h ^ ^ 

Js = (CT ■ u) 

2 2m*i ^ ' 



d__d_ 
dx dx 



The partial derivatives act on the expressions to their right/left respectively (indicated 
by the arrows). For the scattering state (A. 3) we then obtain for the spin current (A.5) 
inside lead i 



(n'cr')Gi 

where {na G j, j = i). For the scattering state (A. 4) we find the corresponding 
expression {na e j, j ^ i) 



I&nai^ e i) 



2m* 



E '^'la',n.(^)f 



(n'(T')ei 

Since every channel is populated according to the Fermi-Dirac distribution f{E; ji^ of 
the respective contact g, the total spin current in lead i then reads 

N 



h{x e i) 



m 



dE 



(A.6) 



-1 r 
— I dE 
47r 



f{E;i,,)B^'{E) + Y,mi,,)n'''{E) 



where 



<t'=± 
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and 

K:l'(E) = Y.J2\tZn'AE)\\ (A.7) 

KAE) - J2J2Kn'AE)\" . (A.8) 
Since S{E) is unitary, the following relation holds true: 

{n'a')ei {n"a")eg 

Then it is straightforward to show that 

mE)+j2n''iE)^o, 

In view of Eq. (A. 6) we eventually find for the spin current in lead i 

h{xei)= — j dEY,[f{E;fi.)-f{E-f,,)]T^'\E). (A.9) 

Although equilibrium spin currents can locally exist in systems with SO interactions as 
shown for 2DEGs [46] and in mesoscopic systems [47], Equation (A.9) clearly shows that 
in thermal equilibrium (/Xj = n) the spin current inside leads without SO interactions 
and magnetic fields vanishes. This absence of equilibrium spin currents in the leads 
has been shown for systems with preserved time-reversal symmetry [48]. However, here 
we show that it is more generally valid for any coherent conductor, since we did not 
make any assumptions about symmetries of the scattering region in the course of the 
derivation. 

We note that the relation (5) derived in Ref. [49] that allows for equihbrium spin 
currents, has to be regarded as incorrect, arising from an improper treatment of the 
back-reflection, missing the term 

- J^dE /(E;/x,)2 {R^^^^E) - R^.^E)) (A.IO) 

in comparison with Eq. (A.9). This issue has been addressed by Nikolic et al. [50]. 
However, their argumentation that the term (A.IO) has only to be included in 
equilibrium, i.e. for energies up to the lowest chemical potential of the N terminals, 
seems questionable. If the spin current is evaluated in a lead connected to a reservoir 
k with Ilk > IJ'i {iJ'i being the lowest chemical potential of any of the N reservoirs), the 
full expression (A.9) has to be used. Furthermore, a simphfied version of Eq. (A.9), 



4:71 , 

where transport at zero temperature and energy-independent transmission probabilities 
were considered, has been used in recent publications on the mesoscopic spin Hall 
effect [51]. 
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Appendix B. Derivation of symmetry relations for spin dependent 
Landauer transport at finite bias 

Here we derive symmetry relations for a two terminal setup as used in this Paper. We 
generalize related expressions from Ref. [52] to finite bias and arbitrary spin quantization 
axis. 

If the Hamiltonian of the total system of scattering region and leads is invariant 
under certain symmetry operations V, we can relate the elements of the scattering matrix 
even in different rocking situations. As an example we derive the symmetry relations 
stemming from the symmetry operator V = —iCRxRu<^z, where C is the operator of 
complex conjugation, R^. inverts the x-coordinate, Ru changes the sign of the applied 
voltage (±^7o ^ =F^o) and cr^ is the Pauli spin operator. 

Generalizing Eq. (A.l) of Appendix A the eigenfunctions inside the leads are given 

by 

with the spin eigenstates 

defined with respect to the quantization axis 

^ sin 6 cos (j) ^ 
u — sin 9 sin 
y cos 9 y 

The effect of V on the eigenstates (B.l) is to change an incoming state in the left lead 
(L) in one rocking situation into an outgoing state of the right lead (R) in the other 
rocking situation and vice versa. Furthermore the position of the spin on the Bloch 
sphere is changed from {9, (p) into {9, — + 7r), and the amplitude of the state is complex 
conjugated. 

On the other hand, the action of V cannot change the scattering-matrix, since the 
Hamiltonian is invariant under the action of V. Therefore the following relation holds 
true: 

n'e(LUR) o-'=±l 

Here, the mode index n is related to mode n of the opposite lead by means of the 
symmetry transformation. Comparing Eq. (B.2) with the inverse of Eq. (A. 2), 

n'G(LUR) cr'=±l 

and using the unitarity of the scattering matrix, S^^ = = (S*)*, we find the symmetry 
relation for the scattering amplitudes 

sI^J^'^He, tUo) = S^:,f.{E, ±U,). (B.3) 
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Similar relations can be obtained for any other symmetry operator that commutes with 
the Hamiltonian H. For the other two operators used in this paper, Rr^RuUz and 
RxRyRuCz, the above procedure can be applied accordingly. It yields 

Si^tZ^ (E, TUo) = 5^41. {E, ±Uo) (B.4) 

for [H, RxRuO'z] — and 

StTJ^E, TUo) = PnPn'Si^SAE, ±Uo) (B.5) 

for [H, R^RyRijaz] — 0, where p„ = (—1)""^ is the parity of the eigenfunction Xniv) in 
Eq. (B.l). 
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